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\ Abstract. It is well known through the work of Majumdar, Papapetrou, Hartle, 

and Hawking that the coupled Einstein and Maxwell equations admit a static mul- 
tiple blackhole solution representing a balanced equilibrium state of finitely many 
point charges. This is a result of the exact cancellation of gravitational attraction 
and electric repulsion under an explicit condition on the mass and charge ratio. 
The resulting system of particles, known as an extremely charged dust, gives rise 
to examples of spacetimes with naked singularities. In this paper, we consider the 
continuous limit of the Majumdar-Papapetrou-Hartle-Hawking solution modeling 
a space occupied by an extended distribution of extremely charged dust. We show 
that for a given smooth distribution of matter of finite ADM mass there is a con- 
\ tinuous family of smooth solutions realizing asymptotically flat space metrics. 
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The purpose of this paper is to estabhsh the existence of an infinite family of smooth 
O \ solutions to the statically coupled Einstein and Maxwell equations within the Majum- 

dar [20] and Papapetrou |24j metric describing the gravitational and electromagnetic 
fields generated from an extremely charged distribution of cosmological dust. 

Consider a system of massive particles carrying electric charges of the same sign. 
^ \ In Newtonian theory, the inter-particle gravitational attraction and Coulomb repul- 

se ' sion both follow inverse square-power laws so that a perfect cancellation of the forces 

is reached when the ratios of masses and charges of the particles satisfy a balanc- 
ing condition. Under such a condition, a system of particles, as a distribution of 
charged dust, is referred to as extremely charged. Thus, due to the balanced forces, 
a static star of any geometric shape may be formed with an extremely charged dust. 
It is obviously important to know whether the same conclusion may be reached in 
the context of general relativity when the Einstein equations are coupled with the 
Maxwell equations. The earliest study of the Einstein-Maxwell equations was carried 
out by Reissner and Nordstrom [35] who obtained a static solution to the coupled field 
equations in empty space, which corresponds to the gravitational field of a charged, 
non-rotating, spherically symmetric massive body and generalizes the Schwarzschild 
blackhole metric. Since the Reissner-Nordstrom solution concerns only an isolated 
single particle, it gives no clue to the above extremely charged dust question. At 
roughly the same time as Reissner and Nordstrom, Weyl |36j classified all Einstein 
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metrics induced by any static axially symmetric distribution of matter and charge. 
The work of Weyl work did not generate as much interest as those by Schwarzschild, 
Reissner, and Nordstrom, due to its geometric restriction of axial symmetry. A 
dramatic development came in 1947 when Majumdar [20] and Papapetrou [24J in- 
dependently extended Weyl's work to include an arbitrarily distributed system of 
extremely charged particles. In particular, the Majumdar-Papapetrou formulation 
gives rise to the electrostatic solutions of the Einstein-Maxwell equations without 
any symmetry restriction, which marvelously generalizes the Newtonian theory to 
the relativistic situation of multiple charged blackholes [11]. It has also been argued 
that, although astrophysical bodies are electrically neutral within a reasonable ap- 
proximation, the Majumdar-Papapetrou solutions may provide simple quasi-static 
analogies for complex dynamical processes preventing asymmetries in gravitational 
collapses or collision of blackholes [T^] . More recently, it has been recognized that the 
Majumdar-Papapetrou solution allows a natural extension to all higher dimensions 
[T8| [19] which is known to be important to issues concerning unification of gravity 
with other forces in nature [211 [22ll26ll27l[28] . 

Within the Majumdar-Papapetrou framework of an extremely charged discrete 
dust model, the coupled Einstein-Maxwell equations reduce to a Poisson equation in 
the vacuum space so that its solution represents a multiply distributed point charges 
with masses [TT]. When a continuous distribution of extremely charged matter is 
considered, the field equations reduce to a semilinear elliptic equation. 

However, in the literature for the continuous case, much attention has been given to 
dust models with spherical star or shell structures [3l|9l[l0l[T2l[T3l[Tl[l6l[I71[l9l[3l] 
with the focus on explicit construction of solutions. There has not yet been any 
study of the existence problem of the governing elliptic equation when the matter 
distribution is arbitrarily prescribed nor of the transition process from the continuous 
model to the discrete model. The purpose of the present paper is to establish an 
existence theory for this important gravitational system. 

In the next section, we follow [12] to lay out the governing elliptic equation for 
extremely charged dusts with a continuous matter distribution, within the framework 
of Majumdar-Papapetrou. In Section 3, we consider the condition to be imposed on 
the mass or charge density which ensures a finite ADM mass. In Section 4, we 
use the methods in Ni [23] to establish an existence theory for the equation. We 
show that, in fact, the Majumdar-Papapetrou equation allows a continuous monotone 
family of finite ADM mass solutions labeled by their asymptotic values. In Section 
5, we derive some asymptotic estimates for the solutions obtained which lead to the 
expected asymptotic flatness of the gravitational metrics. In Section 6, we develop an 
existence theory using an energy method [5, 4j that requires much weaker conditions 
on the decay rate of the mass or charge density function. In Section 7, we present 
some nonexistence results. In particular, we show that when the asymptotic value 
of a conformal metric is larger than some Po > 0, there can be no solution (whereas 
there are solutions when the asymptotic value is less than /3o). In the special situation 
when the mass or charge density function is radially symmetric, the set of permissible 
positive asymptotic values is precisely of the form {0,/3o]. This result complements 
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the conditions in our existence theorems presented in Sections 4 and 6. Finally in 
Section 8, we draw our conclusions. 

2. The Majumdar-Papapetrou metric 

Consider the d = n + 1 dimensional Minkowski spacetime of signature ( — !-■■■+). 
Following \TE\, the metric element is given by 

(2.1) ds^ = 5'^ydx'^dx'^ = —Vdt^ + hijdx'^dx^ ^,1^ = 0,1, ■ ■ ■ ,n, i, j = 1, ■ ■ ■ ,n. 

Here and in the sequel, since we are interested in static solutions, all field variables 
depend on the space coordinates (x*) only. The universal gravitational constant and 
the speed of light are both taken to be unity. Using to denote the covariant 
derivative, the electromagnetic gauge potential, and F^j, = V^Ai, — Vj^A^ the 
electromagnetic field, we can write the electromagnetic energy-momentum tensor as 

(2.2) AttE^, = F^F,^ - ^g^.F^pF^^". 

Moreover, let be the underlying relativistic matter fiuid (i-velocity and pe and 
Pm the associated electric charge and mass densities, respectively. Then the current 
density is expressed as 

(2.3) = PeU^, 

and the material energy-momentum tensor is written as 

(2.4) T^, = p^U^U, + 

where M^^, is the perfect-fiuid, pressure-p, stress tensor given by 

(2.5) M^,=p{U^U, + g^,). 

With the above notation, along with the Einstein tensor Gfj,^ = R^i, — \g^uR rep- 
resented in terms of the Ricci tensor R^^, and Ricci scalar curvature R, the coupled 
Einstein and Maxwell equations are 

(2.6) G^, = M{T^, + E^,), 

(2.7) VyF^"" = AtxJ^. 

As usual, the Weyl-Majumdar-Papapetrou ansatz in higher dimensions is generalized 
to be 

(2.8) A, = <P51, U, = -VV5l. 

Thus, the conservation law V^(T^'^ + E^'^) = 0, which is the relativistic analogue to 
the Euler equation, reads [HI IB] 

(2.9) dip + ^{p„, + p)diV - ^pedi^ = 0, i = l,---,n. 



which takes, after substituting the Weyl-type relation V = V{(j)), the compressed 
form 

(2.10) ^,p+(^{pm+p)^-^y^^ = 0, 2 = l,---,n. 
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On the other hand, the higher dimensional generahzation of the Majumdar-Papapetrou 
condition obtained in HHI reads 



(2.11) V^(0)= (a± 



n 



for some positive constant a and the extremely charged perfect fluid condition [18 
which follows from the compatibility of the Einstein and Maxwell equations, is 



yV/n-2 n 
2.12 Pe = 4-— -pm + -P 

V \n — 1 n — 1 

As a consequence of (12.111) and (12.121) . one arrives [18j at 



2{n-2)( n 
(2.13) + 

In view of ^LW) . (IZTTj) . and (I2A3|) . it is seen that there holds (n - 2)Vd,p = pdiV, 
which, upon an integration, gives us the relation [18] 



(2.14) p = CV^, 

where C is a constant. In a perfect-fluid star, we have p = at the boundary surface 
of the star. Thus, we arrive at C = 0. In other words, p = everywhere and the 
fluid becomes a dust. In particular, we are led to the following general dimensional 
extremely charged dust condition [TS] extending that of the Majumdar-Papapetrou 
solution given as 



2(n-2) 

2.15 Pe = ±\- T^P-^- 

V n — 1 

In the bottom situation with = 3, we have the classical condition for the extremely 
charged cosmological dust 

(2.16) \pe\ = Pm- 

Choose U > and set 
(2-17) V=^^. 

Following [21] and [H], we rewrite the metric element (12. ip diagonally as 

(2.18) ds^ = -fJ-Mt^ + U^dijdx'dx^. 

Then, under the condition (I2.15p . the Einstein-Maxwell equations are reduced to the 
following single nonlinear elliptic equation [18] 

(2.19) AU + 87r(^^-^^Pm{x)U^^ =0, x G M". 

A solution of this equation will give rise to gravitational and electric forces gener- 
ated from the extremely charged cosmological dust distribution. In the subsequent 
sections, we aim at constructing solutions of the equation. 
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3. An ADM mass calculation 

Recall that the ADM mass [U [29] of the space (M", {hij}) housed within the 
Minkowski spacetime of the metric element fl2.ip is given by the formula 

(3.1) "iADM = 7T7 \-\ 1™ / {dihij - djhii)u^ dSr, 

where dSj. is the area element of dB^ {B^. = {x E M" | |x| < r}), v denotes the 
outnormal vector to dBr-, ook is the surface area of the standard unit sphere S^, and 
the metric {hij} is asymptotically Euclidean satisfying 

(3.2) = c5,, + 0(|a;r("-2)), d^g,, = 0(|x|-(""i)), dAg,, = OHxD, 

for I a; I — i- oo with c > 0. Inserting (I2.18P into (13.11) . integrating, and using (I2.19p . we 
obtain 



1 /" / 4 2(ra — 4)^^2{3-nl \ 

3.3 = / —t^^U—^ + > ^ U—VU^ dx. 

^n^i jRr. \{n - 1) (n-2)2 J 

From (13. 3p . we clearly see that, when n > 4, ^aDM — ^ "^ADM = if 
only if = and the conformal metric is flat, U =constant. This statement may 
be viewed as a much simplified (Majumdar-Papapetrou) version of the positive mass 
theorem [251 EOl ISH |32l |33l [37] . Curiously, though, (13. 3 p is not as clear in the classical 
dimension n = 3. Thus, we need to take a different path to recognize "^aDM- 
Since (13. 2 p is translated into 

(3.4) f/ = f/oo + 0(|x|-("-2)), d,U = 0(|x|-("-i)), d,d,U = 0{\x\-^), 

as |x| — )■ oo for some Uoo > 0, we have, in view of (I2.19P again, and the decay estimates 
stated in (13.40 . 

1 f 4-ndU 

m\T\\/\ = — : lim / f/"-2 -— dS*,. 



lim 



{n — 2)w„_i r^oo 

4 — n 



dBr 



4-n 

- u^-' 



dU r iizudU 



(3.5) = ^ USo' I Pmi^)U^^ix)dx, 

{n-l)UJn-l Jun 

which is indeed positive-definite and vanishes only when = and the metric 
becomes globally flat. Thus, we have arrived at a rather transparent (simplified) 
version of the positive mass theorem within the Majumdar-Papapetrou metric in all 
dimensions, n > 3. 

In view of (13.41) and (13. 5p . we obtain the finite integral requirement 

(3.6) / Pm{x) dx < oo. 
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In the subsequent sections, we will construct solutions of the Majumdar-Papapetrou 
equation fl2.19p . observing the condition ()3.6I) . 

4. Existence of solutions by a sub-supersolution method 

Following the above discussion, we rewrite the Majumdar-Papapetrou equation 
over M" (n > 3) in the compressed form 

(4.1) Au + K{x)u^ =0, xgM", 

where K{x) is a smooth nonnegative function in M" which is not identically zero and 
lies in L^(M"). That is, 

(4.2) / K{x) dx < oo. 
So it is natural to assume that 

(4.3) K{x) = 0(|x|-^) as \x\ oo 
where i is large enough. For existence, a convenient condition is 

(4.4) i>2. 

From 04.31) and (14.41) . we may assume that there is a function Kq^v) > (r = \x\) 
such that 

(4.5) K{x) < Koir), r > 0; Ko{r) = 0(r-^) as r ^ oo; £ > 2. 

We may adapt the method of sub- and supersolutions developed in [23]. For this 
purpose, consider 

(4.6) Urr + - -Ur = -KqMu^"^, r > 0, 

r 

(4.7) m(0) = a > 0, Urifd) = 0. 

It is well known that (14. 6p - (14. 71) can be solved locally over [0, R) where R > is 
finite or infinite. Integrating (14. 6p - (14. 71) . we get 

(4.8) r'^'^rir) = - [ p""'^ Ko{p)u^ (p) dp, < r < R. 

Jo 

Hence u{r) decreases in [0,R) if u{r) remains positive-valued. 
Following [23], we consider a test function v defined by 

(4.9) ^W = ^YT^' 
Then v{0) = 1, Vr{0) = 0, and 

(n-1) f 2(a + l)r2\ 2a 

4.10 Vrr + - ivr = -[n- ) ' r-— , f > 0. 

Thus, when 

n — 2 

(4.11) n-2(a+l)>0 or a < — 

V is superharmonic, Av < 0. We shall maintain the condition (14. lip in the sequel. 
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Let u be the unique local solution of fl4.6p - fl4.7p which remains positive over its 
maximal interval [0,R). That is, the solution exists and stays positive in [0,R) for 
some i? > so that u{R) = if i? is finite. 

Define ip{r) = u{r) — av{r), r G [0, R). Then ip{0) = 0, ^9.^(0) = 0, and ip satisfies 

(n-1) , ^ ^ ( 2(a + l)r2\ 2aa 



(l + r2) J(l + r'^)a+i 
2(a + l)r^\ 2aa 



a 



l+r2) y(l + r2)«+i 
2{a+iy\ 2a , , ^ 



(4.12) = g(r), r <R. 

We want to make Q{r) > for all r > 0. To do so, we set 

(4.13) 2(a + l)<£ 

for ^ > 2 given in (14. 5p . Combining (14. lip and (I4.13p . we may assume 

(4.14) a < imin{n-2,£-2}. 

In (14. 5p . we may assume without loss of generality that there is a constant > 
such that 

(4.15) K,{r)<^^^^ r>0. 

In view of (I4.1ip -( l4.15p . we can find some ao > such that 

(4.16) <5(r) > 0, r>0, ae(0,ao)- 
From (I4.12P and (I4.16p . we have as before, 

(4.17) r"~Vr(^) > / P"~^<5(p) dp > 0, < r < i?. 

In particular, <^(r) increases for < r < i? and (/'(O) = implies 

(4.18) lim (y9(r) > or lim u{r) = lim if{r) + a lim v{r) > 0. 

This result establishes R = oo. In other words, there is some ao > such that for each 
a G (0, ao) the problem (14. 6p - (14. 71) has a unique solution, say Ua{r) over < r < oo 
so that fair) = Ua{r) — av{r) increases and Ua{r) stays positive but decreases for all 
r > 0. Consequently, we have the positive finite limit 

(4.19) lim Ua{r) = lim Lp{r) = (3a> 0. 



r—^oo 



As a solution of Au + KqU"-'^ = in R", Ua is a supersolution of (14. ip because the 
property Ua > and (14. 5p imply that 

n n 

(4.20) - Aua = Ko{\x\)uS-' > K{x)ua-\ x eW. 

Set Va{x) = Pa- We see that is a subsolution of (14. ip since K{x) > 0. However, 
we have seen that Va < Ua- As a consequence, we have 
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Lemma 4.1. Equation ( [^. j| j has a solution u satisfying 

(4.21) lia = Vc,<u<Uc,<a, x G M" 

(4.22) u{x) -> (5a as \x\ oo. 



Moreover, u is the unique minimal positive solution of in the sense that if w is 
any other solution satisfying w(x) — )■ /3q, as |x| —> oo, then w > u. 

Proof. We briefly sketch the proof since it is standard. We first solve 

n 

-Auj+i = KuJ~^ in Bii{0) 
Uj+i = /3a on OBr^O) 

Uo = Va 

and observe that 

Va < Uj < Vj^i < Ua (respectively w). 

(Here we used that w > Pa', see section 7 for a simple proof.) Hence ur = Mmj^ooUj 
solves 

-Aur = Kul-' in 5^(0) 
ur = /3a on dBR{0) 

Finally, in the limit as R tends to infinity, we obtain the unique minimal solution u, 
of fll.ip satisfying 

Va ^ u < Ua (respectively w). 

□ 

As a by-product of 04.191] and 04.211) . we have the immediate result 

(4.23) u{x) — )■ /3a as |x| — > oo. 
Next observe that 

(4.24) = (^-pj Ua, < /3 < 
satisfies 

(4.25) ul^^>P, ul^\x) ^ (3 as \x\ oo, 
and 

2 

(4.26) - A«i« = (^j-^ Ko(«f )^ > i^(wf )^. 

Combining ( I4.24p -f l4.26p . we see that we have obtained a sub- and super solution 

pair, consisting of v^a'^ = /3 and u^a \ satisfying v^a^ < Ua^ for each < P < In 
particular, using the above procedure, we see that there is a minimal positive solution 
u^^\x) for any < /3 < (3a which depends continuously on /3. 

In summary, we have 
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Theorem 4.2. Suppose that the function K{x) > satisfies 

(4.27) K{x) = 0{\x\~^), as \x\ ^ oo, £ > 2. 

Then there is a positive number /3o > such that the Majumdar-Papapetrou equation 
Iji4-1\ ) has a continuum of minimal positive solutions {u^^^ | < /3 < /3o} satisfying 

(4.28) lim m(^)(x) = f3, m^^) > (3. 

|a;|— >oo 

Moreover u^^'^ depends continuously in [5 and is monotone increasing in (5. 
In Section 7 we will show that (5q < oo. 

We next consider the asymptotic decay estimates of the solutions obtained in The- 
orem |1]2] as — > oo. 

5. Asymptotic estimates 

We need the following somewhat standard result concerning a Newton-type poten- 
tial integral. 

Lemma 5.1. Define the function v{x) by the Newton-type integral 
(5.1) v{x)= [ r^^dy, 



\y — x\'^ 

where f{x) = 0(|a;|^''') as \x\ — ?■ oo for some constants 7 and a < n and 



(5.2) / \f{x)\dx<oo. 
Then we have 

(5.3) v{x) = 0(|x|^'') as \x\ — )■ 00, 

where S = min{cr, 7-l-cr — n}, provided that the exponents 7 and a satisfy the condition 

(5.4) 'y + a > n. 

Proof. For fixed x G M", x 7^ 0, we decompose into = f2i U U (^3 where 

fli = \ye 



nn 



\X\ 

Y 



^2 = <yeW 



y — x\ < 

\x\ 

— ^ \y — x\ < 2\x\ 



^3 = <ye 



nn 



\y — x\ > 2\x\ 



which results in v{x) = vi{x) + V2{x) + v^lx) with 
(5.5) v,{x)= [ r^^dy, 2 = 1,2,3. 



In,, \y - a^r 

On the other hand, using the condition on f{x), we have the bound 

(5.6) s s tttIpf - rTFr- 
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in view of |x| — \y\ < \x\/2 for y G Qi. Therefore, 

f dy 



\vi{x)\ < 



(5.7) = 0(|x|-(^+"^")) as \x\ oo. 

For 2/ G 1^2, we have 



\v2{x)\ < ^ I \fiy)\dy 

^\x\<\y~x\<.2\x\ 



\X\ 



< ^ [ \f{y)\dy 

(5.8) = 0(|x|-'') as|a;|^oo. 

Furthermore, using ||x — y| — |x|| < \y\, we have 



\V3{X)\ < C 

< c 

< C 





y-x\>2\x\ |y-a;|'^(l+ bl^) 



\y-x\>2\x\ \y - a;|'"(l + [|y -x\ - |a;|]T) 
dy 



,-x\>2\x\ |Z/-X|-(1+ |y-x|V2^) 

(5.9) = 0(|xr(^+"-")) as |x| ^ oo, 

using I ly — a;| — |x| I > \y — x\/2, which is the same as stated in fl5.7p provided that 
7 + cr > n holds. 

Hence, summarizing fl5.7l) - fl5.9p . we see that (15.31) is estabhshed. □ 

We now turn our attention to (14.11) . Let m be a positive solution of (14. ip produced 
in Theorem 14. 2[ Then there is a constant /3 > such that 

(5.10) lim u{x) =13. 

\x\—^oo 

Let w{x) be given by the Newton potential associated with (14. ip . i.e., 

n 

5.11 w{x) = — 7W?r ~\ 

n{n - 2)iln Jru \x - yY" ^ 

where f2„ = 27r"/^/nr(n/2) is the volume of the unit ball in M". Applying Lemma 
15.11 with 0" = n — 2 to w{x) and noting (I5.10p . we see that, with £ = 7, the conditions 
(14. 4p and (15. 4 p coincide. Therefore, we obtain 

(5.12) w{x) = 0{\x\~"') as |x| 00, m = min{n - 2, £ - 2}. 

On the other hand, since h = u — w is a. bounded harmonic function in view of 
(I5.10p and (15.120 . it must be a constant, which must be given by (15.100 . In other 
words, we have established the relation 
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Differentiating fl5.13p . we liave 

(5.14) = 0(|a;r("+^)) as |a;| ^ oo, i = l,---,n, 

in view of Lemma 15.11 again. 

In order to get tlie decay estimate for D'^u near infinity, we may be tempted to 
differentiate fl5.13p twice and apply Lemma 15.11 again. However, sucli an approacli 
will not work because of the technical restriction a < n for the Newton-type integral 
(15.1 p stated in Lemma I^TTl To overcome this difficulty, we set \x\ = R, assume that 
R is sufficiently large, and rewrite (I5.13P as 

(5.15) u{x) = (3 + U{x) + V{x), 
where 



(5.16) Uix)= Tix-y)fiy)dy, V{x) = V{x - y)f{y) dy, 

J\y\<2R J\y\>'iR 

with T{z) = l/n{n - 2)1]„|2|2-" and f{y) = -K{y)u^{y). Then \didjV{z)\ < 
First, we have by using \x — y\> \y\ — R and fiy) = 0{\y\^^) for \y\ — )■ oo. 



mV{x)\ < / \d^d,T{x~y)f{y)\dy 

J\y\>2R 

< c j 

J\y\>2R F " Vr 



< TC j |?/|-(^+")di/ 

J\v\>2R 



l?/l>2i? 

(5.17) = CiR-\ 2,j = I,-- - ,n. 

In order to estimate the second derivatives of U{x) in (I5.15p . we extend /(y) such 
that f{y) = for \y\ > 2R and use Lemma 4.2 in |7j to represent didjU{x) as 

d,d,U{x) = [ did,T{x-y){f{y)-f{x))dy-f{x)[ d,T{x ~ y)u,{y) dSy 

J\y\<3R J\y\=^R 

(5.18) = f/i(x) + f/2(x), 2,j = l,-- - ,n, 

where {vjiy)) is the outnormal vector at y over the sphere {\y\ = 3R} and dSy is the 
area element. 

In view of |a; — y| > \\y\ — |a;|| = 2R for \y\ = 3R, we have 

\U2ix)\ < \f{x)\ [ \d,T{x-y)u,{y)\dSy 

J\y\=3R 



< C|/(x)|i?-("-i) / dS, 

J\y\=3R 

(5.19) < CiR'^. 



12 JOEL SPRUCK AND YISONG YANG 

To estimate Ui{x), we let flk {k = 1, 2, 3) be defined earlier and write 

Uii^) = E /" d,d,T{x-y){fiy)-f{x))dy 
k=i Jnkn{\y\<m} 

(5.20) = Wi{x) + W2{x) + W-i{x). 

We impose the additional natural decay condition on K as follows, 

(5.21) diK{x) = 0{\x\-^^^^^) as|x|^oo, i = l, 
In view of fl5.2ip . f l5.14p . and the relation / = Ku^^ , we have 

(5.22) dif{x) = 0{\x\'^^+^^) as|x|->oo, i = l,---,n. 

Since for y G Qi, we have \y\ < \y — x\ + \x\ < 3R/2 < 2R, we see that Wi{x) satisfies 
the estimate 



\W,{x)\ < / \d,d,T{x-y){f{x)-fiy))\dy 

J {\x-y\<R/2}n{\y\<2R} 
J\x-y\<R/2 F y\ 
J\x-y\<R/2 1-^ y\ 

where ^ = x + t{y — x) for some t G [0, 1]. Hence |^| > |x| — ||/ — x| > R — R/2 = R/2 
and |V/(^)| < Csi?"^^"*"^^ for some constant C3 > 0. Inserting this result into fl5.2c 
and integrating, we obtain 

(5.24) \Wi{x)\ <C^R-^. 
To estimate the decay rate of W2{x), we first note that 

(5.25) / ld^d,Tix-y)fiy)\dy<-^ [ \fiy)\dy. 

J §<\y-x\<2R iin-ft J^n 

Moreover, in view of f{x) = 0{R^^), we have 

2 



\d,d,Tix-y)fix)\dy < 0(i?-^)p^ / dy 

§<\y~x\<2R 



^§<\y-x\<2R' ^ \ /\ \ 'Q^Ri^Jr 

(5.26) < 2^''0{R 



Combining (I5.25P and fl5.26p . we obtain 

(5.27) \W2{x)\<C5R-"^'''^^'''^. 

For Wslx), we make a similar decomposition. First we have the estimate 

(5.28) / lO^O,T{x-y)f{y)\dy<-^ [ \fiy)\dy. 

J\y~x\>2R ^ J^n 
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Furthermore, by virtue of f{x) = 0{R^^) again, we find 

/ \d,d,Tix-y)fix)\dy < j 

(5.29) < (-\ 0{R-'). 



Combining fl5.28p and (15.29p . we arrive at 

(5.30) |Vr3(a;)| < C75i?-™"^^'">. 

Consequently, we may now collect f l5.19p . f l5.24p . f l5.27p . and f l5.30p and recall the 
decompositions fl5.18p and fl5.20p to obtain the decay estimate 

(5.31) \didjU{x)\ < Celxp"'^^^^'"^ as |x| ^ oo, t,j = !,■■■ ,n. 

Finally, inserting fl5.3ip and fl5.17p into fl5.15p . we see that the second derivatives 
of the solution u{x) satisfy the same decay estimates as U{x). 
Summarizing the above discussion, we can state 



Theorem 5.2. Let u be a solution obtained in Theorem with u{x) l3 > as 
\x\ — )■ oo when the charge density function K satisfies the conditions and {4-4)- 
Then u enjoys the decay estimates 

(5.32) u{x)-l3=0{\x\-'^), Du(a;) = 0(|a;|~(™+^)), \x\^oo, 

where m = min{n — 2,i — 2} and D denotes a generic derivative. Moreover, if 
the function K satisfies the additional condition h5.21\) . then D'^u obeys {D^u){x) = 

0(|x|-("^+2)) ^ 

The asymptotic flatness condition (13. 4p may be ensured with the assumption £ = n 
such that m = n — 2. 

6. Existence of solutions by an energy method 

In this section, we shall develop an existence theory for the Majumdar-Papapetrou 
equation (14. ip using an energy method. 

We assume that K{x) satisfies (14. 3 p with i > 1. If the constant /3 > defined 
below is positive, we assume in addition that (14. 2 p is fulfilled. 

Choose a sequence {9^} of positive constants decreasing to zero and define the 
Hilbert space 

nl = {v e H\W) : [ {\Vv\^ + 9kV^)dx<oo}, 
with inner product 

{v,w)k= / (Vv ■ Vw + 9kvw) dx . 



and norm \\v\\k = (f,f)^. For later use, define also the Hilbert space Hk as the 
completion of the smooth functions of compact support under the inner product 



{v,w)k 



/ (Vf ■ Vw + Kvw) dx . 
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We will be interested in the subset 

Ak = {v enl:v>0, ||i;||fc = l} . 

Fix vo e C°°(M") nnln L°°(M"), ||vo||o < l,vo > 0. We define a sequence Vk G 
Ak, k > 1 inductively as follows: for a fixed constant (3 > 0, let Wk+i G 'Hl_^_l be the 
unique solution of 

(6.1) -Aw + ek+iw = K{vk + P)^ 
and set Vk+i = \k+iWk+i where Afc+i = p;^^^- Then 

(6.2) Vk+i E Ak+i and - Avk+i + 9k+iVk+i = \k+iK{vk + /3)^ . 

To see the legitimacy of the above procedure, we inductively assume that Vk G "H^ 
and Wfc ^ some k>l. 

Recall the Gagliardo-Nirenberg inequality 

n — 2 

(6.3) ( / v^dx] " <C \Vv\'^dx. 

As a consequence of (14. 2 p and (16.31) . we see that K{vk + P)"^ G L^(R") since 

2 r 2n 2 

K\vk + f3)^dx < Ci K\v;:-' + (3^) dx 



n-2 



< C2 / \WvkVdx] +C3 / is:dx<oo. 

Using the Lax-Milgram theorem, we see that fl6.ip has a unique solution, say Wk+i, 
in 'H\j^i, as claimed. 

We need to show that Wk+i > 0. To see this, we consider the weak form of (I6.ip . 
that is, 



(6.4) 



/ {ywk+i-V(l) + 9k+iWk+i(t))dx= I K{vk + P)"-^(t)dx, (f) e Ul^^. 
Jr" Jr" 



Let m > and observe that = (— m — Wk+i)^ G 'Hl_^_-^. Inserting this choice of 
into (16.41) gives 



6k+iWk+i(l)dx > / K{vk + 13)^^-2 (j)dx >0, 

{«)fc+i<-m} Jr" 

which is impossible unless the set {wk+i < —m} is of measure zero. Since m is 
arbitrary, we have shown that Wk+i > a.e., which establishes the legitimacy of the 
construction of the sequence {vk} as stated. 
We define 

$(1;) = / K{v+ + /3)^dx . 
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Note that for v G Ak-, we have the uniform bound 



2n-2 2n — 2 , 



^{v) < 2"-2 / K{v'^-2 + /3 "-2 ) dx 

< clp"-^ [ Kdx+\\K\\Ln^RnJ [ \Vv\'dx) 
(6.5) < c(f3'^ [ Kdx + \\K\\Ln(^^^y 



with a uniform constant C > 0. 

Lemma 6.1. ^{v^) < ^{vk+i), /c = 0, 1, . . . 

Proof. From f l6.ip . we have 



(Vffe+i ■ V(wfc+i - ^fc) + 6'fc+it;fc+i(t;fc+i - ffc)) 



= Afe+i / K{vk+i + P)ivk + /3)"-2 dx - <l>(t;fc) 

72 — 2 

(6.6) =^--^A,+i($(^;,+i)-$(t;,)). 
However, the left-hand side of (16. 6p satisfies 

(Vffc+i ■ V(wfe+i - Vk) + Ok+iVk+i{vk+i - Vk)) dx 

(6.7) > l{\\v,+i\\U, - \\vk\\l + {Ok - Ok+i)\\vk\\l2^^r.^) > 0. 



Corollary 6.2. Assume that (5 lies in the interval 
(6.8) < /3 < /3o = 



□ 



Then there is a positive constant C such that the sequence {A^} lies in the interval 
C\\KlJ§^ (||i^||Li(R") + \\K\\LHmY^' < < Ao = Hvor' (l - £ 
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Proof. From (16 ■2p and Lemma [6. ![ we have 



1 > / {Vvk+i ■ Vvk + 9k+iVk+iVk) dx 
= Xk+i(Hvk)-f3 f K{vk + /3)^^dx 



2n-2 



> Afc+i|$K)-/3 / Kdx] ^Vk)—^ 

> \k+Mvo)(l-4-]. 



/3i 



Similarly, in view of Vk+i G Ak+i, the inequality fl6.3p . and the bound (I6.5p . we also 
have 



1 = Xk+i / Kvk+i{vk + /3)"~^ dx 



< 



Tl-2 

2n 



n-2 



< \k + lC\\K\\lZ,^„s (||i^||il(Rn) + ||-ftr||i;,n(]Ku)) 



2 



for a uniform constant C > 0. □ 

In order to prevent the sequence {Xk} from diverging to infinity or trivializing at 
zero, from now on, we assume /3 lies in the interval (16. 8p as stated in Corollary 16.21 

Corollary 6.3. || bfe+i - H+i < ^Xk+i{'^{vk+i) - '^{vk))). 
Proof. We have 

^ll^^fe+l -^^fcllfe+l = {Vk+l,Vk+l-Vk)k+l + ^{\\Vk\\l+i-\\Vk+l\\l+i) 
< {Vk+l,Vk+l - Vk)k+1 

n — 2 

by I^M) and Corollary [O since ||t^fc||fc+i < \\vk\\l = 1 = □ 

We next show the sequence {vk} is uniformly bounded and obtain some (non- 
optimal) decay estimates. 

Proposition 6.4. There is a uniform constant C such that for k > 1, 

(^^)vk{y)<C{l + \y\y- where a = mm , - e) , (0 < £ < if /3 > 0, 

^>,(l/)<C(l + |i/|)-^ z//3 = 0. 
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Proof. Fix a ball = Bji{y) (of radius R > and centered at y E M") in M". Then, 
by standard local elliptic estimates, we have 

(6.9) supt;fe+i < C (R~'^\\vk+i\\Lp{Ba) + R^~'^\\K{vk + /3)^^\\LoiBii)) , 

Br ^ ' 

"2" 

n 

P>0, ? > 2" 

We first observe that ffc(x) — > as |x| — ?■ oo for all A; > 1. In fact, since K 
is bounded and satisfies (14.21) . we see that K{yQ + /3)^^ G Lp(]R") for any p > 1. 
For p > n, elliptic theory applied to (16. ip implies Vi G iy^'^(]R"). Thus Vi{x) — )■ 
as |x| — )■ oo. Iterating this argument, we see that Vk+i{x) — > as |x| — )■ oo for 
A; = 0,1,2,---. 

To prove (i), let i? = 1, p = and choose ?/ such that supf^+i = Vk+i{y)- Then 
applying (16. 3p and Vk E Ak {k > 1), we obtain 

(6.10) supvk+i<c(^l + e(^j {vkY^^dxy +C{e) 

for ^ < q < n (for example q = Therefore, we find from (I6.10p and (16. 3p the 
bound 

supffc+i < C ^ + e{supvk) (^j Wfc""' dx^ +C{e 

(6.11) < C{l + 6supvk + C{e)). 
Now fix e = 2^ to obtain 

(6.12) sup f fe+i < C + ^ sup Wfc 

for a uniform constant C. Iterating, we find 

(6.13) sup Vk+i <2C + 2"(^+^) sup Vq, 
proving (i). 

To prove (ii), let 2R = \y\, p = ■^^2- Then, from (16. 9p . part (i), and the boundedness 
of $(va:), we have 

v.^M < c[\y\-^"-^^ + \yr^.-'-^^^) 

(6.14) = cf|y|-("f') + |yr('-^+"^)^ 



If / > n, we let g = oo and the result follows immediately. If n > £ > ^, then 
/ - 2 + ^ > ^ for q sufficiently close to | and the result follows. If 1 < / < |, 
choose ^ = f — -^11 completing the proof of (ii). 

Finally if /3 = 0, we observe that the decay estimates of part ii. are still valid and 
if we have shown that Vk{y) < C(l + |?/|)~", a < then by the argument of part 
(ii), we find 



..^,(,)<c(|yr(^) + i,r('-+^+"-^)) 
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SO that Vk+i < C(l + \y\)~'' with b = + - ^) > Hence in a finite 

number of steps we arrive at the decay rate for a (different) uniform constant 
C. □ 

Remark 6.5. The dependence of the constant C on vq comes from the sup norm 
estimate of part (i). In particular, from fl6.13p . we see that this dependence will 
disappear in the limit as k tends to infinity. Also if we only assume \ Vvo\'^dx < 1 
(instead of HfoHo < 1), then, from (16. 6p and (16. 7p with /c = 0, we have 



2n-2 



1„ „ ,,2 



(6-15) > -2^o||^^o|lL2(M")- 

Hence, in view of ( 16.15^ and Lemma [6.11 we have 

(6.16) Hvk+i) > Hvi) > <^{vo) - Ceo\\vo\\l2^^ny 
This will be used in the proof of Theorem 16.71 below. 

We are now ready to state and prove our main existence theorem of this section. 

Theorem 6.6. The sequence {vk} has a subsequence converging locally uniformly in 
to a smooth positive function v such that u = v + P is a finite energy solution 

of 

(6.17) -Au = XKu^ %nW, 

(6.18) u{x) (3 as \x\ — )■ oo, 

(6.19) / \Vu\^dx < 1, 



where A > 0. If /3 = 0, u = 0{\x\ ^)) while if (3 > 0, u - (3 = 0(|x|~") where 
a = min (^^, — e), {0 < e < ^"~^'* /). Moreover, all estimates are independent 

of sup Vq by Remark lUTSi 

Proof. By Lemma W/U and (16. Sp . {$(ffc)} converges to a finite limit, say M. Since the 
sequence {vk} is uniformly bounded in Tix H L°°(]R"), we can choose subsequences 
{ffc-}, {ffci+i} which converge weakly in nL°°(]R"). By Corollary 16.31 both subse- 
quences converge weakly to the same limit, say v. Using that each v^+i satisfies (16.21) 
and Afc+i stays uniformly bounded from above and below (Corollary 16.20 . Proposi- 
tion [6]1] and standard elliptic regularity estimates give that the convergence is locally 
uniformly in C2+"(M") and that u = v + ^ e C°°(M") satisfies (l6TrD - (l639D and 



2n-2 



(6.20) / K{x)u^{x)dx = M, 

where A > is a limiting point of the sequence {Aa,}. The decay properties oiu = v+13 
follow from Proposition 16.41 if /3 > and using a bootstrap argument as in part (iii) 
of Proposition 16.41 in case /3 = 0. □ 
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We next use Theorem 16.61 to solve the global constrained variational problem asso- 
ciated with equations f l6.17l) . fl6.18p . Consider the variational problem 

(6.21) Mp := max < $(f ) | v 

As shown earlier, Mg is always well-defined and the question is whether Mg is 
achieved. 

Theorem 6.7. Let K be a smooth nonnegative function satisfying K = 0{\x\~^), i > 
1. If (3 > 0, assume Iji4-^ addition. For (i in the interval Ii6. 8\) . the variational 
problem Mg has a finite energy solution v satisfying ( [g. 1 7p , ^.18\) . ( fg. with the 
decay properties described in Theorem 1 6'. 61 

Proof. Let v^"^^ be a maximizing sequence for Afg, that is ^{v^"^^) > — m = 
1,2,---. Without loss of generality we may assume v^""^ E C°°(]R") n H^iR"") fl 
L°°(M'^), v'^"^^ > 0, $(l;*^'")) > ^ . Given m > 1, we choose in our iterative scheme 

^0 = ("^ll^''™'' |Il2(M"))~"^5 '^o'"'* — construct the sequence {v^^^} converging 

to the solution v^"^^ (say) of Theorem 16.61 Then from Lemma 16.11 (see f l6.16p and 
Remark 16. 5p , 

<|.(t;(™)) > $(t;i™))>$(l;(™))-C0o||t^^'"^|li2(K„) 
C + 1 

(6.22) > Mg . 

m 

Since \Vv^"^^\'^dx < 1, the sequence {v^"^^ is also a maximizing sequence of solu- 
tions satisfying 

-Aw^'™) = A('")ir(t;(") + P)^ . 
By our previous arguments, a subsequence converges locally uniformly in C^+"(M") 
to a nonnegative solution v of —Av = \K{v + /3)"-2 for some A > which is a 
maximizer of the variational problem f l6.2ip . □ 

A solution w of the original Majumdar-Papapetrou equation (14. ip may then be 
constructed from a solution u of the nonlinear eigenvalue problem fl6.17p by setting 



n-2 



(6.23) w = \—u. 



7. Nonexistence results 
For the Majumdar-Papapetrou type boundary value problem 

(7.1) Au + K{x)u^ = 0, M>0, xgM", 

(7.2) lim u{x) = /3, /3 > 0, 

\x\—^oo 

where K{x) > is a smooth function which is not identically zero and satisfies 

(7.3) [ K{x)dx<oo; K{x) = 0{\x\-^), i > 2 or £ > 1, 
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we have shown that a solution exists when /3 is not too large. Here we note that there 
may be no solution when /3 is sufficiently large. Without loss of generality, we choose 
the origin so that 

(7.4) K(0) = msixK(x) > 
and set 

(7.5) Ko{r) = min{i^(x) \ \x\=r,x e M"}, r > 0. 

Lemma 7.1. If the boundary value problem consisting of ( [7. ip - fT!^ has a solution, 
then the modified problem, 

(7.6) Av + Ko{r)v'^ = 0, Mm v{x) = /3, v>0, 

X'|— >-oo 

has a radially symmetric solution. 

Proof. Let m be a solution of f l7.ip - fl7.2l) . We ffist show that u > (3. In fact, since 
u{x) — )■ /3 as |x| — )■ oo, we see that w = max{0,7 — m} = (7 — u)+ is of compact 
support for any < 7 < /3. Multiplying (17. ip by w integrating by parts, we have 

(7.7) - / |VMpdx= / K{x)u^^wdx>0, 

which proves that the set {u < /?} must be empty as claimed. 

On the other hand, using -K'o(^) = -K'o(la^l) < K{x),x G MJ^ {\x\ = r), we have 

(7.8) -AuyKou"^. 

In other words, f + = m is a supersolution of (17. 6p . Taking v~ = (3, we see that v~ 
is a subsolution of (17. 6p and v~ < v~^. We can now construct a radially symmetric 
solution V , of (17. 6p satisfying t>~ < t> < v~^. This solution is just the minimal positive 
solution. As before, we solve 

-Avj+i = Kov;^' in Br{0) 
Vj+i = /3 on dBii{0) 
Vo = v~ 

and observe that since the Laplace operator commutes with rotations, the sequencejfj} 
is radially symmetric by uniqueness. Moreover, 

Hence VR^r) = limj^oo Vj{r) solves 

n 

-AvR = Kov^-' in Br{0) 
VR = /3 on 95^(0) 

Finally, in the limit as R tends to infinity, we obtain a radially symmetric solution, 
V, of (17. 6p satisfying t>~ < t> < t>+. □ 

Proposition 7.2. Let v{r) be a radially symmetric solution of the boundary value 
problem ( [7. 6[ ) where Kq{0) > 0. Then v is apriori bounded by a constant M (indepen- 
dent of (3). 
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Proof. Integrating (17. 6p over the interval (0, r), we see that v = v{r) decreases. Hence 

(7.9) v{0) = max{v{r) | r > 0} 

If the proposition is false, then we may find a sequence of radially symmetric solutions 
of (17. 6p . say {vj}, satisfying 

(7.10) Mj = Vj{0) ^ oo J = 1, 2, ■ ■ ■ . 
Now for fixed j = 1, 2, ■ ■ ■ , define the rescaled function 



Then Q < Wj <1 = Wj (0) and 



(7.12) - A^,- = M. Kvj = M. M3:)v;-' (x) = iro(x)w;-' . 

where A denotes the Laplace operator with respect to the variable x. 

Using elliptic theory and passing to a subsequence if necessary, we see that the 
sequence {wj} converges (locally uniformly in C^"*"") to a function w on with 
w(0) = 1, which satisfies 

(7.13) Am; + Ko(0)w^ = 0. 

On the other hand, however, by Theorem 1.1 in [6j, we have w = 0, which gives us a 
contradiction. □ 

Corollary 7.3. Suppose that the function K satisfies ( [y.^[ ). // the constant /3 > 
in \7.2\} is sufficiently large, the boundary value problem consisting of ( [7. i[ ) and Ij?.^ 
will have no solution. 

Proof. If (I7.ip - (l7.2p has a solution, then (17. 6p also has a radially symmetric solution 
as described in Lemma 17.11 Hence by Proposition 17. 2[ /3 < M. □ 

Theorem 7.4. There exists Pq > so that there is a (minimal) solution of the 
boundary value problem ( [7. i| ) and ^1.2) for any < /3 < /3o and no solution for 
P > Pq. If K[x) is radially symmetric about the origin with -ft'(O) > 0, there is also a 
solution for Pq. 

Proof. Define /3o be the supremum of the asymptotic values /3 > such that the 
boundary value problem (17. ip and (17. 2p has a solution. As we have shown in the proof 
of Theorem 14. 2[ there exists a minimal positive solution u^^\x) for < /3 < /3o. Now 
assume that K{x) = Ko{\x\) is radially symmetric about the origin with Kq{0) > 0. 
Then by Proposition [721 the minimal positive solutions u^^\r) with asymptotic value 
(3 are all uniformly bounded independent of /3. Hence ^'^^^(r) = lim/3_^^(, M*^'')(r) is a 
solution of (17. 6p with asymptotic value /3o- □ 

Next we consider a "shell" -star type solution of the equation (17. ip . To motivate 
our discussion, recall that (17. ip allows a many-particle solution of the form [TT] 

N 

(7.14) u{x) = c + Y. ^ 



n- 
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where c, /xi,--- are positive constants. Such a solution determines the metric 
induced from a system of extremely charged particles of masses /xi, ■ ■ ■ , /xat located 
at the points pi, ■ • • ,pAr G M" so that the mass density is expressed as a sum of the 
Dirac distributions concentrated at pi, - ■ ■ ,Pn- In particular, u{x) — )■ oo as x — )■ 
{j = l,---,N). 

Similar, an idealized "shell-star" type solution [10] represents a situation in which 

(7.15) K{x) = Ko{x)6{F{x)), 

where S = {x E M" | F{x) = 0} is the shell which can be realized as the boundary 
surface of a bounded domain Q, i.e., S = dfl, and Ko{x) > 0. Suggested by the afore- 
discussed many-particle solutions, we may imagine that f2 = fii U . . . U fi^ is the finite 
union of simply connected compact domains and look for a solution of fl7.ip - fl7.2l) in 
M" \ satisfying the exterior boundary condition 

(7.16) m(x) ->cx), x-^S = dn. 
However, one easily sees that this is hopeless. In fact. 

Lemma 7.5. There is no function u that is superharmonic (Au < 0) and bounded 
below inW^ \Q and satisfies ^7.10^ . 

Proof. Choose the origin to lie in the interior of VL. By adding a large constant to 
M, we may assume m > 7 > in M" \ Vt. Then by the (weak) maximum principle, 
u > 1^1^-2 in M" \ f2 for any positive C > 0. This is clearly impossible. □ 

8. Conclusions 

In this paper, we have studied the Majumdar-Papapetrou equation describing an 
arbitrary continuously distributed extremely charged cosmological dust in static equi- 
librium governed by the coupled Einstein and Maxwell equations in a general (n+ 1)- 
dimensional spacetime and established several existence theorems for the solutions 
of the equation assuming that the ADM mass is finite. Among the results are the 
following. 

(i) Assume that the mass or charge density vanishes at infinity at a rate no slower 
than 0(|a;|~^) for some positive number £. When £ > 2, a positive solutions ap- 
proaching a positive asymptotic value at infinity can be constructed by a sub- and 
supersolution method. When the condition is relaxed to £ > 1, a positive solution 
with the same properties can be obtained by an energy method. The obtained solution 
and the structure of the equation allow us to establish the existence of a continuous 
monotone family of solutions labeled by their positive asymptotic values at infinity 
which are necessary for achieving asymptotic flatness. 

(ii) Specific asymptotic rates of a solution obtained above have been established 
under various decay assumptions on the mass or charge density. In particular, if 
i = n and the first derivatives of the mass or charge density decays like 0(|x|~*-""'"^^), 
then the conformal metric of a solution obeys all the decay rate assumptions in the 
classical definition of asymptotic flatness at infinity. 

(iii) The set of possible asymptotic values of the conformal metrics determined by 
the Majumdar-Papapetrou equation is a bounded interval. In other words, there is 
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a positive number /3q > 0, such that there is no solution u satisfying the asymptotic 
condition u{x) — )■ /3 as |a;| — )■ cxd when P > Pq and there is a solution when P < Pq. 

(iv) There is no such "shell-star" solution of the type that the solution gives rise to 
an asymptotic fiat space and approaches infinity near the shell which is the boundary 
of a bounded domain so that the mass density is concentrated at the shell given by a 
Dirac distribution. 
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